Tracing the ancetsry of galaxies 
Proceedings IAU Symposium No. 277, 2010 
Carignan, C, Freeman, K. C. & Combes, F. eds. 



© 2010 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



The impact of ISM turbulence, clustered 
star formation and feedback on galaxy mass 
assembly through cold flows and mergers 

Leila C. Powell 1 , Frederic Bournaud 1 , Damien Chapon 1 , Julien 
Devriendt 2 , Adrianne Slyz 2 and Romain Teyssier 1 

1 Service d'astrophysique, CEA, Orme des Merisiers, Gif-sur-Yvette Cedex, France 

email: leila.powell@cea.fr 
2 Oxford astrophysics, Denys Wilkinson Building, Keble Road, 0X1 3RH 



Abstract. 

Two of the dominant channels for galaxy mass assembly are cold flows (cold gas supplied via 
the filaments of the cosmic web) and mergers. How these processes combine in a cosmological 
setting, at both low and high redshift, to produce the whole zoo of galaxies we observe is largely 
unknown. Indeed there is still much to understand about the detailed physics of each process in 
isolation. While these formation channels have been studied using hydrodynamical simulations, 
here we study their impact on gas properties and star formation (SF) with some of the first 
simulations that capture the multiphase, cloudy nature of the interstellar medium (ISM), by 
virtue of their high spatial resolution (and corresponding low temperature threshold). In this 
regime, we examine the competition between cold flows and a supernovae(SNe)-driven outflow 
in a very high-redshift galaxy (z ~ 9) and study the evolution of equal-mass galaxy mergers 
at low and high redshift, focusing on the induced SF. We find that SNe-driven outflows cannot 
reduce the cold accretion at z ~ 9 and that SF is actually enhanced due to the ensuing metal 
enrichment. We demonstrate how several recent observational results on galaxy populations 
(e.g. enhanced HCN/CO ratios in ULIRGs, a separate Kennicutt Schmidt (KS) sequence for 
starbursts and the population of compact early type galaxies (ETGs) at high redshift) can be 
explained with mechanisms captured in galaxy merger simulations, provided that the multiphase 
nature of the ISM is resolved. 

Keywords, methods: numerical, galaxies: evolution, galaxies: formation, galaxies: ISM, galax- 
ies: star clusters 



1. The impact of feedback from individually- resolved SNe on a 
protogalaxy with a multiphase ISM 

The Nut simulations (presented in detail in Powell et al ( 201Clh ) probe the evolution 
of a protogalaxy embedded in the cosmic web at z ps 9, at very high resolution (0.5pc 
physical in the densest regions) using the resimulation technique with the AMR code 
RAMSES (|Tevssierll2002l) .The main simulat ion includes cooling , SF, S Ne feedback (kinetic 
feedback modelled as Sedov blastwaves dDubois fc Tevssierl [20081) ). metal enrichment 
and a UV background (jHaardt h, Madaulll996T ). A galactic wind develops, filling a large 
fraction of the volume out to ~ 6r v i r with hot gas, which is punctured by cold filaments 
of the cosmic web that propagate all the way to the disc. SF rates are high for the galaxy 
halo mass (1M Q /yr for a 10 9 M Q halo) and we report the formation of star-forming clumps 
with typical masses of 10 6 M Q . In order to isolate the impact of the SNe-driven galactic 
outflow, we undertake a parallel analysis of a control run with no SNe feedback. 
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Figure 1. Right: The global 
average SF efficiency of the main 
galaxy versus redshift in the 
feedback (solid line) and no-feed- 
back (dashed line) runs. Left: 
Ratio of the net inflow in the 
clumpy component to the net in- 
flow in the filamentary compo- 
nent versus redshift for the feed- 
back (solid line) and no-feedback 
(dashed line) runs. 



1.1. Cold inflows versus hot outflows 

We establish temperature and density thresholds that separate out the filaments (0.1 ^ 
p ^ 10 atoms cm~ 3 , T < 2 x 10 4 K), the clumps (i.e. condensed gas, such as the main 
disc and satellite galaxies) (p > 10 cm~ 3 , T < 2 x 10 4 K), cold diffuse gas (p < 0.1 cm~ 3 , 
T < 2 x 10 4 K), warm diffuse gas (p < 10 cur 3 , 2 x 10 4 K < T < 2 x 10 5 K) and hot diffuse 
gas (p < 10 cm -3 , T > 2 x 10 5 K). We are then able to measure and compare the inflow 
and outflow rates in each component; we calculate the mass flux through spherical shells, 
centred on the halo centr e, out to the virial radius. The full results of the inflow/outflow 
analysis are described in iPowell et al (l2010h . but here we recall the main conclusions: 



Inflows The most striking result is that the far-reaching wind does not impede the ac- 
cretion of gas onto the central galaxy. We find that the filaments are the dominant supply 
of gas to the disc, and since the gas in the filaments is dense and highly supersonic, it 
is not surprising that the wind has not been able to halt its progress. This suggests that 
for galaxies with host haloes of this mass (10 9 M Q ) that are predominantly fed by cold 
streams, SNe-driven winds are unlikely to be able to limit accretion and therefore reduce 
SF in these systems, as is commonly suggested. 

Outflows Note that we associate the galactic wind with the warm and hot diffuse com- 
ponents, defined above. We find that while the wind is far-reaching (extending to around 
6r v ; r by z = 9), the mass outflow rate is only about 10 per cent o f the inflow ra te; we do 
not find a massive wind as seen in low redshift observations (e.g. lMartinlll999l ). We also 
report some outflow in the cold, diffuse component, which we attribute to gas that was 
in the region around the filaments being entrained by the wind and swept out; an effect 
that is seen in observed low-redshift winds. 

1.2. The resulting star formation 

Surprisingly, we find the star formation rate (SFR) in the Nut simulation is a factor of 
« 10 higher than that in the control run (without SNe) at z = 9. We attribute this to 
the metal enrichment which results from the inclusion of SNe feedback. We find that the 
global average efficiency of SF is very high in both the Nut simulation and the control 
run. We use an efficien cy of 1 per cent to control SF on parsec-scales in RAMSES, as found 
in the observations of iKrumhofz and Tanl (|2007l ). However the global efficiency is ~ 10 
per cent as shown in Fig. [1] (left). We find that this efficiency is correlated with merger 
activity; in Fig. [1] (right) we show the ratio of inflow in the clumpy component to that 
in the filamentary component. High values of this quantity indicate that mergers are 
supplying a significant fraction of the gas being accreted; it is clear that these instances 
coincide with peaks in the global SF efficiency of the main galaxy (left). 

These characteristics of the SF can be explained as follows. As evidenced by the star- 
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Figure 2. Right: Gas density 
pdf before (solid) and during 
(dashed) the merger and at the 
peak of the starburst (dot dash). 
Left: The KS relation for discs 
(solid ) and starbursts (dashed) 
from iDaddi et all (|2010h ■ and se- 
lected points for simulated galax- 
ies before the merger (diamonds), 
during the starburst (triangles) 
and after the merger (squares). In 
the latter stage the galaxy has 
been morphologic ally quenched 
(|Martig et aUl2009D . 



forming clumps seen in t he Nut simulation a t z = 9, and supported by clumps captured 
in ot her simulations (e.g. Agertz et al. 20091) and observations of clump-cluster galaxies 



(e.g. lElmegreen fc Elmegreenl 120051 ). all at z = 3, it seems probable that SF occurs in 



dense clumps at high redshift. We postulate that the correlation of the global SF efficiency 
with the relative importance of merger activity reflects the fact that mergers can in duce 
fragmentation of the ISM, resulting in a more clumpy medium ( Tevssier et al. 2010l ) and 
more gas becoming eligible for SF. We examine the process of merger-induced SF in 
detail in the next section. 

2. Merger-induced star formation at high resolution in a cloudy ISM 

Merge r-induced clustered SF has been demonstrated in simulations of the Antennae 
system ( Tevssier et al. 2010l ). We are undertaking a similar investigation, but with a 
whole sample of mergers with different orbital parameters, in order to assess the relative 
importance of this type of merger-induced SF. We have simulated 5 1:1 mergers, with 
AMR code RAMSES with a maximum spatial resolution of 5pc, allowing us to resolve 
densities up to 10 6 cm~ 3 . The sample is currently being extended both with additional 
orbits and different mass ratios. Here, we focus primarily on examining the mechanism 
via which SF is induced. 

2.1. Changes in the star formation and ISM structure 

We find that, in some cases, the peak in the SF rate occurs before the final coalesence 
of the two galaxies and the gas density distribution shows an extended, clumpy region, 
rather than a centrally-peaked but otherwise smooth distribution. In order to understand 
what impact the merger has had on the ISM we examine the gas density probability 
density function (pdf) in the left of Fig. [2j Here we can see the evolution of the pdf from 
before the merger (black line), during the initial stages of the merger (blue line) and at 
the peak of the starburst (red line). It is clear that there is an excess of very dense gas 
(p > 10 4 cm -3 ) which increases as we progress through the merger. 

2.2. Insights into observations of merging systems 

It is clear that in order to reproduce the formation of star clusters, which are observed 
both in high redshift clump-cluster galaxies and in lower redshift merging systems, it 
is necessary for hydrodynamical simulations to have a high spatial resolution (and the 
associated low temperature threshold). However, aside from the clustered SF, additional 
effects are apparent in such highly resolved simulations of merging systems which are not 
seen when the ISM is warm and smooth (as in low resolution simulations). 
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• We have captured an excess of dense gas directly in our merger simulations (as 
demonstrated in the left of Fig. HI). It was pro posed (and demonstrated using post- 



processing of simulations) by I Juneau et al (2009) that an excess of dense gas in the ISM 
could be induced by mergers and that, since HCN emission traces denser gas phases than 
CO, this could account for observations of enhanced HCN/CO ratios in ULIRGs. This 
would remove the need to invoke active galactic nuclei activity. 

• The increased fraction of gas mass in this very dense phase (in the form of clumps) 



can al so provide a physical explanation for the two KS relations demonstrated in baddi et al 



( 20101 ). Fig. [2] (right ) shows that as our merger simulations progress, the galaxies move 
up and to the right in the EsFR-Egas plot, reaching the starburst sequence at the peak 
of the SFR, then moving down to below the disc sequence once much of the fuel has 
been consumed. There is some central gas inflow which increases the average gas surface 
density E gas . However, the merging galaxies reach the starburst sequence because Esfr 
is much higher than would be expected from a smooth gas distribution with the same 
value of Egas due to the high efficiency of SF in the very dense gas clumps. 

• Simulations of high redshift mergers (which have a higher gas fraction than those 
discussed above) show that the resolving th e clumpy nature of th e ISM can also have 



a significant impact on the merger remnant (jBournaud et a J l2010h . There are two main 



effects when the ISM is multiphase: 1) the merger remnant is much more compact and 2) 
there is a much smaller surviving/reformed disc component. The first conclusion provides 
a possible formation mechanism for the population of compact ETGs at high redshift, 
while the latter casts doubt on the likelihood of discs surviving/reforming after mergers 
without some external source of gas (e.g. cold accretion). 

3. Conclusions 

We demonstrate that many insights into galaxy mass assembly via both cold flows and 
mergers can be gleaned from simulations with sufficient resolution to capture the cloudy 
nature of the ISM and clustered SF. We find that SNe-driven outflows cannot reduce the 
cold accretion onto dwarf-mass galaxies at z m 9 and that the SFR is actually higher 
due to the ensuing metal enrichment. We describe how high-resolution galaxy merger 
simulations enable us to explain several observational results (e.g. enhanced HCN/CO 
ratios in ULIRGs, a KS sequence for starbursts and high-redshift compact ETGs) as they 
reveal physical processes not captured previously. 
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